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Abstract 

By applying the particle-number projection to the finite-temperature BCS theory, the 5-shaped 
heat capacity, which has recently been claimed to be a fingerprint of the superfluid-to-normal 
phase transition in nuclei, is reexamined. It is found that the particle-number (or number-parity) 
projection gives 5-shapes in the heat capacity of nuclei which look qualitatively similar to the 
observed ones. These 5-shapes are accounted for as effects of the particle-number conservation on 
the quasiparticle excitations, and occur even when we keep the superfluidity at all temperatures 
by assuming a constant gap in the BCS theory. The present study illustrates significance of the 
conservation laws in studying phase transitions of finite systems. 

PACS numbers: 21.90.-tf,21.60.-n,05.70.-a 
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While phase transitions are striking phenomena in inhnite systems, it has been difficult 
to establish them in hnite systems such as atomic nuclei, because their signatures are often 
obscured by the quantum fluctuations [l|. Since most nuclei have the superfluidity in their 
ground states, i.e. at zero temperature, the superfluid-to-normal transition for increasing 
temperature has been discussed theoretically Recently, high precision measure¬ 


ment of nuclear level densities has been implemented p, in which the level densities are 
extracted from the 7 -ray data with the help of the Brink-Axel hypothesis. Converting the 
microcanonical information to the canonical one, they have found that R-shapes appear in 
the graphs of the heat capacity C* as a function of temperature T. It has been argued that 
the R-shapes are a hngerprint of the superfluid-to-normal phase transition, since such S- 
shapes occur if continuity is taken between heat capacity in the superfluid phase (described 
by the constant-T model) and that in the normal fluid phase (described by the backshifted 
Bethe formula). Based on this picture, they estimated critical temperature from their exper¬ 
imental data. Similar R-shapes come out in theoretical calculations including the quantum 
fluctuations 0 , 0 . 

The superfluid and the normal-fluid phases in nuclei are defined within the mean-field 
picture; in the Bardeen-Cooper-Schrieffer (BCS) theory or in the Hartree-Fock-Bogoliubov 
theory. In studying the superfluid-to-normal transition, linkage to the mean-held picture 
should carefully be traced, even when the R-shapes are reproduced by calculations including 
a variety of quantum huctuations. Breakdown of a certain symmetry is often associated with 
phase transitions. In respect to the superhuid-to-normal transition, the particle-number 
conservation is violated in the superhuid phase while preserved in the normal-huid phase, 
within the mean-held theories. However, this is an approximate picture and the particle- 
number conservation is not violated in actual nuclei, restored via the quantum huctuations. 
Since restoration of the symmetry sometimes plays an important role particularly in hnite 
systems, it is desired to investigate how the number conservation ahects the heat capacity 
of nuclei. 

The so-called number-parity forms a subgroup isomorphic to S 2 of the U{1) group ac¬ 
companied by the particle-number. The number-parity projection in the hnite-temperature 
BCS (FT-BCS) theory 0 was developed more than two decades ago 0. The full number 
projection at hnite temperature is much more complicated task. There is a fundamen¬ 
tal difficulty in the variation-after-projection scheme, though an approximate solution has 


recently been suggested 


On the other hand, the number projection in the variation- 

In this paper, we apply 


lumoer prpj 

• 0Hln|. 


before-projection (VBP) scheme was formulated in Refs, 
the particle-number projection as well as the number-parity projection in the FT-BCS the¬ 
ory, and qualitatively investigate ehects of the particle-number conservation in the heat 
capacity of nuclei. By the VBP scheme, we view ehects only of the projected statistics, 
without changing the excitation spectra given by the BCS Hamiltonian, as a step of tracing 
ehects of various quantum huctuations. 

We mainly consider the nuclei. Let us assume the following model Hamiltonian, 
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where k indicates the time-reversal of the single-particle (s.p.) state k. The s.p. state k and 
its energy Sk are determined from the Nilsson model |^, by assuming the quadrupole defor¬ 
mation from the measured E2 strength |^. We take Qp = 22/A and Qn = 27/^4MeV jl^ . 
For the model space, we hrst dehne the Fermi energy ep for each nucleus by the arithmetic 
average between the energy of the highest occupied Nilsson s.p. level and that of the lowest 
unoccupied level, without the residual interaction. We then include all the s.p. levels sat¬ 
isfying \ek — ^fI < 7MeV, for both protons and neutrons. The Qr values are related to the 
model space; with the present choice we reproduce the pairing gaps of 12/-\/]4MeV. 
Although the Hamiltonian in Eq. m is relatively simple, it will be sufficient for qualitative 
study of the pairing phase transition. Note that the critical temperature of the deformed-to- 
spherical shape phase transition is appreciably higher than that of the superfluid-to-normal 
transition, and that the s.p. levels hardly change at T < 1 MeV 0|. 

In the FT-BCS theory, we introduce the auxiliary Hamiltonian 


H' = H'+H/, Hl = Hr-XrNr {T=p,n), 


( 2 ) 


where At stands for the chemical potential and iW the number operator, W = J^ker 
The quasiparticle (q.p.) operators are obtained from the Bogoliubov transformation. 


ttfc = Ukttk 


Vkal, 


with ul + vl = 1, and the q.p. energy is given by 




Qr'^k 
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The density operator in the FT-BCS theory is 

g-ho/T 


Wo = 


Tr(e-^o/T) 


= Ekalak 


( 5 ) 


Here Tr denotes the grand-canonical trace in the model space. The thermal expectation value 
of an observable O is calculated by (O)o = Tr(t()oO)- Equation (jHj) is an approximation of 
H' in the Boltzmann-Gibbs operator by (const. -|- Ho)- The entropy is dehned by 

S = —Tr(ti)o Int&o)- The FT-BCS equation, which determines the Uk and Vk coefficients in 
Eq. (jni), is obtained so as to minimize the grand potential H = {H')o — TS for each T j^. 
The chemical potential A,- is hxed by the particle-number condition {Nt)o = where W 
is the particle number in the model space corresponding to the specific nuclide. 

Introducing the particle-number projector. 


Ft = — / 
27r J —TT 


— TT 
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( 6 ) 




we define the density operator in the projected statistics, 


p p p-Ho/Tp p 

U^p = — -- . (7) 

Tr(FpP„e-^o/T) 

The thermal expectation value in the number-projected statistics is obtained by (0)p = 
Tr(tDpO). The integration over ip in Eq. (jH)) takes account of a certain part of the two-body 
correlations beyond the mean-held approximation. In practical calculations, the ip integral 
is replaced by a discrete sum: 


Pr 


1 
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where M stands for the number of the s.p. states and pm = 27rm/{M -|- 1). If we set M = 1 
instead of the number of the s.p. states in Eq. (jH)), Pt is reduced to the number-parity 
projector [^, 


p' = 1 ^ p-MNr-NP 

^ ^=0,7r 


( 9 ) 


Correspondingly, wp and (0)p become those of the number-parity projection, which we 
shall denote by wp/ and (0)p/ = Tr(t()p/0). In this paper we distinguish the unprojected, 
the number-projected and the number-parity-projected expectation values by the suffices as 
(O)o, {0)p and (0)p/, respectively. 

We here calculate heat capacity by C = d{H)/dT, which is obtained by numerical differ¬ 
entiation of (H) for various T, in practice. In Fig. ^ the heat capacities with and without 
projection (Cq, Cp/ and Cp) are depicted for ^®^Dy. There occur two discontinuities in C{T), 
corresponding to the superfluid-to-normal transitions for protons and neutrons. This is con- 
hrmed by Ap and A„, which rapidly vanish at the respective critical temperature Tc. This 
signature to the transition does not disappear by the projection because it is implemented 
after variation. Although such discontinuities are unrealistic, the present purpose is to in¬ 
vestigate effects of the projection mainly at T < Tc. We have Tc ~ 0.5 — 0.6 MeV in the 
present calculation, slightly lower for protons than for neutrons. In the normal fluid phase, 
the projections do not make important differences; difference of the number-parity-projected 
result from the unprojected one is even almost invisible. At T < T^ we hnd that there is a 
certain effect of the projection on C{T), either the number-parity or the number projection; 
an S'-shape appears when we apply the number or the number-parity projection. It is noted 
that, in the number-parity projection, the zero-point of the energy may influence Cp/, since 
(W)p' is displaced from N^. depending on T because of the incomplete projection (after 
variation). To suppress this influence, we use {sk — £f) instead of Ek in the Hamiltonian of 
Eq. dH). 

Since the proton and neutron degrees-of-freedom are separated in the Hamiltonian of 
Eq. (dj), we have C = Cp + Cn, where = d{Hr)/dT. We present Cn(T) for ^®^Dy in 
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FIG. 1: Heat capacity C vs. temperature T for The dashed line is obtained by the FT- 

BCS calculation without projection, while the solid and dotted lines are with the number and the 
number-parity projections, respectively. 

Fig.Efa), so as to simplify discussions. There is no qualitative difference between Cp{T) and 
CniT). It is noticed that, while Cn,o = d{Hn)o/dT increases almost linearly at (0.2 MeV < 
)T < Tc, an F-shape comes out both in = d{Hn)p'/dT and Cn,p = d{Hn)p/dT in 

this temperature region. In order to view effects of the projections more clearly, we also 
show AEn^p = {Hn)p — {Hn)o ^nd AEn^pi = {Hn)p' — {Hn)o as a function of T, in Fig. |2Kb). 
Obviously, the slope in AE^^piT) is equal to the difference of from Cn^, and likewise for 
AEn,p/{T). In this respect, the F-shape in the number-parity-projected case comes from the 
decrease of AEn^p/{T) at 0.2 < T < 0.35 MeV as well as the increase at 0.35 < T < 0.5 MeV. 
This behavior does not change by the full particle-number projection, although AEn,p{T) 
shifts downward to a certain extent, compared with AEn,p'{T). We also plot in Fig. |2Kc) 
the expectation value of the q.p. number (7v„), where In Fig- l2Kc), the 

number-parity-projected result is almost indistinguishable from the fully number-projected 
result at T < 0.5 MeV, and from the unprojected one at T > 0.5MeV. At T 0, we have 
{Mn) ~ 0 as it should be. As T goes up slightly to about 0.2 MeV, excitation to the 1 q.p. 
states gives rise to increase of {H)o, and yields dominant contribution to Cnfi- However, this 
is hctitious, resulting from the violation of the number or the number-parity conservation. 
Since the 1 q.p. states are removed, {H)p' tends to stay at its T = 0 value, giving the 
decrease of AEn,p'{T) and delaying the rise of C'n,p'(T). As T grows further, excitation to 
the 2 q.p. states starts contributing to {H)o, and also to {H)p'. If more and more q.p.’s are 
excited, the q.p. number distributes broadly, and the lack of the odd number-parity states 
is no longer important. Then AF^„p/ becomes vanishing again, and Cn,p' also approaches 
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Cn,Q- This feature is inherited in the number-projected result. Thus the occurrence of the 
S'-shape in C'(T) is accounted for as effects of the number (or number-parity) conservation 
on the q.p. excitations. 

In the proton degrees-of-freedom are almost the same as in ^®^Dy. Cn{T), AEn{T) 

and (A/”„) in ^®^Dy are presented in Fig. El In the graph of C„(T), we hnd an S'-shape 
somewhat similar to that in ^®^Dy, as a result of the projection. However, the S-shape 
in the projected result is less conspicuous, particularly around T = Tc(~ 0.6MeV). It is 
remarked that similar even-odd difference is observed in the experiments and in more 
realistic calculations SQ- Moreover, we view small but non-vanishing values at low T 
(0.05 < T < 0.2MeV) in the projected results (see inset to Fig. Eta)). These even-odd 
differences are again accounted for in terms of the q.p. excitation picture. Since the neutron 
number is odd in ^®^Dy, the q.p. number should be 1 at T 0, which is not taken into 
account in the unprojected result. Therefore, AE^ is higher by about A„(Ri 1 MeV) at 
T 0 than in ^®^Dy. As T goes up, AE^^p' decreases, because excitation to the 2 q.p. states 
predominantly contributes to (.H)o, which is eliminated in {H)pi. For T > 0.5 MeV, many 
q.p. states mix up and the lack of even number q.p. states becomes less important. This 
leads to AEn,p' ~ 0. However, since AEn,p>{T = 0) is higher than in ^®^Dy, AEn^p' ~ 0 
at T > 0.5 MeV implies that AEn,p' does not go up to a great extent at T ^ Tc, making 
the upper bend in the S'-shape of Cn,p'(T) weak. Although AEn,p shifts downward, the 
full number projection gives analogous behavior. The structure in Cn{T) at low T is also 
connected to the q.p. excitation. While there is only a single 0 q.p. state present in the 
even systems, there are several 1 q.p. states having close energy to one another. Hence, 
excitation from the lowest-lying 1 q.p. state to the higher-lying 1 q.p. states is possible even 
at low T. Note that this effect is only the cases for the projected statistics. To illustrate this 
mechanism, the q.p. number corresponding to the three Nilsson s.p. levels adjacent to 
and that for the next nearest six levels (three higher and three lower levels) are separately 
depicted in the inset to Fig. EKc). Although (ACi)p or (ACi)p' stays unity at T < 0.3 MeV, 
we view that excitation among the s.p. states occurs, giving small but non-negligible heat 
capacity. 

As expected from the q.p. excitation picture, the above results are insensitive to nuclide, 
except that C(T) at low T(< 0.3 MeV) for odd nuclei somewhat depends on the s.p. levels 


around £f- This has been conhrmed by calculations for neighboring Dy isotopes [l^- The 
q.p. excitation mechanism is so generic that it should not depend on the nuclear shapes, 
which has also been conhrmed by calculations for Sn isotopes (l .5| . 

The number (or number-parity) conservation thus gives rise to the F-shapes in the heat 
capacity and their even-odd differences. Since it is explained within the q.p. excitation 
picture, the S'-shaped heat capacity does not seem straightforwardly linked to the superhuid- 
to-normal transition. To investigate this point further, we try the following calculation: 
instead of solving the FT-BCS equation at each T, we keep using the solution at T = 0 for 
ik and A,-, and therefore for Uk and Vk- By this treatment the nucleus stays in the superhuid 
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FIG. 2: Thermal properties of the neutron part of (a) C„, (b) AEn, (c) {Mn), as functions 

of the temperature. See Fig. ^ for conventions. 
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(a) 



T [MeV] 


FIG. 3: Thermal properties of the neutron part of (a) Cn, (b) AEn, (c) (A/’n). Cn at 

0.07 < T < 0.35 MeV is amplified in the inset to (a). In the inset to (c), (Z^fceSPin (red 

lines) and (X]fcGSP 2 „ (blue lines) are also presented as well as the total q.p. number (black 

lines), where SPl^ is composed of the Nilsson s.p. level with ep and the neighboring two (one 
higher and one lower) levels, while SP2„ consists of the next neighboring six (three higher and 
three lower) levels. See Fig. ^ for other conventions. 





















T [MeV] 


FIG. 4: Heat capacities for ^®^Dy: (a) Cn and (b) C, with and A,- fixed to be the T = 0 values. 


phase at any T. The T-dependence enters only via the explicit one in the Boltzmann-Gibbs 
operator. Resultant Cn and C = Cp + Cn are shown in Fig. HI for ^®^Dy, and in Fig. El for 
^®^Dy. We have the S-shapes even though the superfluidity is kept at any T. The even-odd 
differences in C are also viewed in this artihcial model. 

The observed S'-shapes in C(T) are deduced from the energy-dependence of the level den¬ 
sities j^. We have conhrmed that the present calculations with projection are qualitatively 
consistent with the measured level densities, whether or not Hq is T-dependent. However, 
the present results indicate a mechanism of the S'-shapes in C'(T) different from the picture 
assumed in Ref. [^, in which the S'-shapes emerge via continuity between the two phases. 
Still, one should not immediately conclude that the S'-shapes shown in the present work are 
irrelevant to the transition. The value in the present FT-BCS calculations is close to the 
temperature region where the S'-shapes appear, as in the experiments, which may suggest 
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FIG. 5: Heat capacities for (a) Cn and (b) C, with and A,- fixed to be the T = 0 values. 


that the S'-shapes might have correlation to the transition, even though their connection 
is not straightforward. For instance, suppose that the mixing of many q.p. states, which 
causes the upper bend of the S'-shape in C'(T), also drives the phase transition in the BCS 
approximation, the S'-shapes are indirectly correlated to the transition. It has been shown 
that the S-shapes takyjlace (and the discontinuity at critical temperature is washed out) in 
realistic calculations |3|,|5|, in which various correlations are taken into account. However, 
in practice, relation of the S-shapes in C{T) to the superfluid-to-normal phase transition 
has not yet been clear enough. It is desired to inspect carefully how individual quantum 
fluctuations affect the heat capacities of nuclei. We just point out here that the number 
conservation should play a certain role in producing the S-shapes, even though it may not 
be a full explanation. It is emphasized that effects of the conservation laws should not be 
discarded in discussing phase transitions in hnite systems. 
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In summary, we have reexamined the heat-capacity of nuclei, applying the particle- 
number projection to the hnite-temperature BCS theory. Since the projection is carried 
out after the variation, the sharp discontinuity at critical temperature remains which is not 
observed in experiments. While the S-shapes in the heat capacity have been claimed to be a 
hngerprint of the superfluid-to-normal phase transition, it is found that the particle-number 
projection gives S'-shapes in the heat capacity of nuclei analogous to the observed ones, apart 
from the discontinuity. The even-odd difference in the heat capacity is also produced by the 
projection. Except low T part of odd nuclei, the number-parity projection gives similar heat 
capacity to the full number projection, if the s.p. energies are appropriately shifted. These 
S'-shapes are accounted for in terms of the quasiparticle excitations, and occur even when 
we keep the superfluidity at all temperatures by assuming a constant gap in the BCS theory. 
Although the observed S-shapes could still correlate to the phase transition, their relation 
should be inspected carefully. The even-odd difference is also understood in the context of 
the quasiparticle excitations, in which the particle-number (or the number-parity) conser¬ 
vation plays a crucial role. The present study reveals signihcant role of the particle-number 
conservation in the heat capacity of nuclei. 
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